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ABSTRACT 
The  use  of  Laplace  Transforms  has  long  been  one  of  the  primary 
methods  utilized  in  solving  linear  differential  equations  but  recently  the 
extension  and  application  of  this  method  to  the  solution  of  basic  non- 
linear systems  has  been  proposed.    This  latter  technique  is  presented 
and  expanded  to  include  various  non-linear  functions  .    The  applicability 
of  this  method  to  systems  of  varied  form  and  complexity  is  then  explored 
and  the  validity  of  the  derived  solution  investigated. 
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I.     INTRODUCTION 

Laplace  Transforms  have  been  typically  applicable  and  limited  to 
linear  systems.    However,  Baycura  [Ref.   1  and  2]  in  two  recent  papers 
has  proposed  an  extension  of  the  Laplace  Transform  method  to  non-linear 
systems.     By  expanding  the  Laplace  integral  into  an  infinite  series  and 
integrating  by  parts  n  -  times  an  infinite  series  was  developed  for  non- 
linear terms.    Through  a  further  correlation  of  these  derived  series  with 
series  of  linear  transforms  the  transform  for  a  non-linear  term  was  obtain- 
ed.   These  papers  develop  this  technique  for  higher-order  derivatives  and 
powers  of  the  dependent  variables  illustrated  with  several  applicable 
examples  . 

The  purpose  of  this  paper  will  be  to  present  the  general  method  as 
developed  by  Baycura  and  to  extend  this  method  to  varied  and  more  com- 
plicated equations  of  non-linear  systems  .     The  latter  will  be  accomplished 
by  developing  the  transforms  for  the  product  of  derivatives  of  the  depen- 
dent variable  and  the  dependent  variable  raised  to  the  n      power,  and 
derivatives  raised  to  the  n      power.     The  validity  of  the  transforms  will 
then  be  shown  in  several  applicable  examples  .    In  the  examples  a 
computer  method  designed  to  find  the  inverse  Laplace  Transform  will  be 
introduced  and  explained. 


II.     LAPLACE  TRANSFORM  CONCEPT 


A.      THE  METHOD  OF  INTEGRATING  BY  PARTS 

The  following  generally  summarizes  Baycura's  method  of  utilizing 
Laplace  Transforms  on  non-linear  systems  using  the  method  of  integration 
by  parts.     Considering  the  product  of  two  time-dependent  functions,  x(t) 
and  y(t),  which  are  to  be  integrated  using  integration  by  parts  n  -  times, 
obtain  the  general  formula: 


[ xydt  = 


x 


ydt  -  x 


f  ydtdt  +  X-    [J'J*  ydtdtdt  - + 


(_!)n-i  x(n-i)  KJydt11-1^)  dt  +  .  .  . . 


(2-1) 


n  =  1 ,  2  ,  3 ,    ... 
where  x^n-1'  indicates  the  (n   -    1)       derivative  of  x. 

Applying  this  general  formula  to  the  Laplace  Integral,  equation  (2-1) 
can  be  used  to  find  the  series  form  of  the  Laplace  Transform  of  a  function 
x(t)  allowing  y(t)  =  e~s1:.    This  integration  gives  the  series  form  of: 


xe    st  =  x 


00 

.  r 


e"st  dt  -  x  j  J'e_st  dtdt+x     J'J  e_stdtdtdt  -  .  .  .      (2-2) 

9  0°  o 

Evaluating  at  the  upper  limit  it  is  found  that  the  value  of  the  integral 
approaches  zero  as  t  approaches  infinity.    Then  evaluating  at  the  lower 
limit  the  integral  becomes: 


fxe-st   _     x{0}    +  *{0}  +    x(PJ 

jxe  -       g        +      sa    +      ^     +  .... 


(2-3) 


or 


«,r  u\ i  /   \         ^10)  MO)  MO) 

£[x(t)]    =  x(s)    =  -^     +    -V     +    -V    +    •  ... 


(2-4). 


B.      APPLICATION  TO  DERIVATIVES 

Equation  (2-4)  depicts  the  Laplace  Transform  of  x(t)  in  an  infinite 
series  expansion  involving  the  initial  conditions  of  the  function.     Through 
similiar  manipulations  the  transforms  of  the  derivatives  may  be  obtained: 

*  WHS -f *:•.-**    =m     +^     +    ^    +    .... 

(J  o  o  o 

o 

=   sX(s)    -   x(0)  (2-5) 


X[X(t)]  =  fxe-stdt  =  ^   +    =$ 


+ 


=  sy  X(s)  -  sx(0)    -  x(0) 


(2-6) 


C.     APPLICATION  TO  NON-LINEAR  FUNCTIONS 

The  simplest  non-linear  function  to  investigate  is  xa  (t) .     Using  the 
same  technique,  expansion  by  parts  yields: 


£[x2(t)]    =  x 


xe   st  dt  -  x 


I 


xe    stdtdt  +  x* 


J'J'xe    stdtdtdt 


(2-7) 


Taking  each  of  the  above  terms  individually  and  again  expanding  by  parts 
obtain: 


x 


Jxe"stdt  =  x(0)  X(s) 


(2-8) 


-x 


xe   st  dtdt 


MO)   X(s) 


(2-9) 


5c 


i-r       -st  ,     ,     ,        x'(O)   X(s' 
JJxe   SIdtdtdt  =     w  5  v  ' 


(2-10) 


The  summation  of  these  terms  then  yields 

x(0)     +     x(0 


X[x3(t)]    =    [x(0)  + 


+ 


X(s)  =  sX3(s) 


(2-11) 


Through  similiar  computations  ,  transforms  of  higher-order  functions 
may  be  obtained.    The  transform  of  the  function  cubed  would  be: 


£[x3(t)]  =  [x(0)    +  ^    + 


mi 

„2 


(2-12) 


Generalizing  it  can  be  stated  that: 


Z  [xn(t)]  =   s11"1  Xn(s) 


(2-13) 


III.     LAPLACE  TRANSFORMS  FOR  NON-LINEAR  FUNCTIONS 


A.      TRANSFORMS  FOR  THE  PRODUCT  OF  A  FUNCTION  AND  ITS  FIRST 
DERIVATIVE 

Transforming  x(t)  x(t)  by  applying  integration  by  parts  the  following 

derivation  is  obtained: 


X[*(t)x(t)]  = 


xxe"     dt  =  x 


xe   stdt  -  x 


xe   st  dtdt 


(3-D 


+  x 


xe   st  dtdtdt  - 


-st 


Let:    Ij.  =  x  I   xe   SIdt 


Knowing  the  transform  of  the  term  under  the  integral,  this  term  then  yields 

Ix   =  x(0)   X(s)  (3-2) 

Similiarly  allow: 


-V     Pl'ir^-St 


lP= -x    J'xe"stdtdt  =  x(0) 


X(s) 


(3-3) 


and; 


la  =  * 


JO 

'Jjxe-Stdtdtdt  =  x(0)  ^ 


(3-4) 


The  summation  of  these  terms  gives: 


I  [xx]   =  Z     In    =    [Ii   +  Ia  +  I3   +    .  ..  3 
n=l 


X 


(s)  [x(0) 


x-(O)         x(0) 

T  T  j  T 


(3-5) 


Simplifying: 


£  [x(t)   x(t)]  =   s2X(s)[-  ^    +^ 


s2     +    s 


2        +    "73       + 


=   s2  X(s)  fx(s)    - 


f]] 


=  sX(s)    [sX(s)    -  x(0)]  (3-6) 

Continuing  in  this  fashion  the  Laplace  Transform  for  x(t)xa  (t)  is 


derived. 


£  [x(t)  x2(t)]  =  Jxx2e"stdt  =  x 


00  OS 

x2e"stdt   -  k' jjx2e"stdtdt 


+  x  f  jj'x3e   stdtdtdt   - 


(3-7) 


Let:  I,    =  x 


x2e   stdt  =  x  |x 


[- 


xe   st  dt   -  x 


-st 


J'xe   SIdtdt  + 


Further  let:   ]1   =  x  [  x  e   stdt  =  x(0)  X(s) 


J2  =   "  x 


r         -stJ+J.        x(Q)   X(s) 
J    xe   bLdtdt  =  — "-1 *-*- 


J3  =  x  ffj"  xe=stdtdtdt  = 


x*(0)   X(s] 


Summing: 


Mo)   .  x(pj 


Ii   =  *  |=1  Jn  =  *«»  x<s)  [x(°)  +     s       + 


f 


and  after  simplifying; 
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Ix   =   sx(0)  X2(s) 
Similiarly  let: 

00 

Is  =  -  x  [f'x2e~stdtdt  = 


(3-8) 


x  [x 


xe    stdtdt  -  x 


Jjxe~stdtdtdt  +  .  .  .1 


Lj_    =  x 


),e-sididt=-^^^ 


L2  =  -  x  jJJ  xe-stdtdtdt  =  -  *MjXlS) 


Summing  again: 

I2=-xf=1Ln   =   -x(0)X(S)[-^-f^ 

=  x(0)  X2(s) 
By  similiar  arguments  obtain: 


xjO] 
„a 


I3  =  x  jjjx2  e"stdtdtdt 


x(0)  X3(s) 

s 


(3-9) 


(3-10) 


Combining  the  above  integral  terms  the  transform  is  then: 

x  tx(t)  x»jwp  -   I.  i„  =  .'  rfw  r-  5121  +  £fil  +  *i°l  +  .  .  . 

=  s2  X2(s)  [s  X(s)  -  x(0)]  (3-11) 

In  similiar  fashion  the  transform  for  x(t)  x3(t)  may  be  obtained  and 
expressed  as : 


-st 


X[*(t)  x3(t)]-     xx3e    aLdt  =  s3X3(s)    [sX(s)-x(0)] 


(3-12) 


Equations  (3-6),    (3-11)  and  (3-12)  provide  the  Laplace  Transform  for 
the  general  form  of  the  function,  x(t)  xn(t). 


11 


£[x(t)   xn(t)J   =  snXn(s)    (£[x(t)]  } 


(3-13) 


B.      TRANSFORMS  FOR  THE  PRODUCT  OF  A  FUNCTION  AND  ITS  SECOND 
DERIVATIVE 

Employing  again  the  method  of  integration  by  parts  to  the  non-linear 

function,  x(t)x(t),  the  transform  may  be  expressed  as: 


03 


03 


X[x(t)x(t)]   =  Jxxe    stdt  =  x  Jxe   stdt  -  x\^xe   stdtdt  +  ...(3-14) 

0  0  0 

Attacking  each  of  the  integrals  separately  as  in  Section  III. A  we  let: 


lx   =  x*  fxe*stdt  =  x-(O)  X(s) 


(3-15) 


Ip  =  -  x 


I,  =  x<«) 


J'  xe"stdtdt  = 


x(0)  X(s] 


]j ■xe-stdtdtdt  '.  xlilSLxM 


(3-16) 


(3-17) 


Therefore: 


(*> 


l  [x<t)x(t)]  =  g=]  ln  =  x(0)  X(s)   ♦  HMJOsl   +  x(4)(o)  X(s)  +  _ 


=   s 


x(s) r. m .  m  +  xio}  +  m  +  m  + ..  i(3-18) 

L      s  s  s  s  s  J 


Collecting  terms  and  simplifying: 

£[x(t)  x(t)]  =   s  X(s)   [s2  X(s)    -   sx(0)    -   x(0)]  (3-19) 

The  derivation  of  the  Laplace  Transform  of  x*(t)x2(t)  may  be  similiarly 

obtained. 


£[*(t)x2(t)]  = 


x  x2e~stdt 


x 


xae~st  dt  -  x 


xae   stdtdt  + 


.  (3-20) 
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Knowing  the  transform  of  x2(t)  again  write  separately  for  each  of  the 
integrals: 


L   =  x" 


rx2e"stdt  =  k*(0)  sX2(s) 


(3-21) 


In     =    "  X 


r  ,,2 


-st 


x2e        dtdt  =  x(0)  X2(s) 


(3-22) 


I,  -  x<*> 


xae-stdtdtdt  =  xl!ki!M 


(3-23) 


Summing  the  above  terms  the  transform  becomes: 

£[x-(t)x2(t)]=E=1In   =  M(0)sX2(s)   +  xV0)X2(s)  +  x(4)(0)  X2(S)  +  ... 


=  s4X3(s)  [-  l 


S  b«+       B       +^+"^    +    -  (3"24) 


J 


X[x(t)  xs(t)>  s2X2(s)   [s2X(s)  -  sx(0)  -  x(0)] 
The  transform  for  k(t)x3(t)  may  be  identically  expressed  as 


(3-25) 


£[x-(t)  x3(t)]  =  jxx3e"stdt  =  kjx3e    st  dt  -  x'fjx3e   stdtdt  +...  (3-2  6) 

0 

o  o 

JC[k'(t)  x3(t)]  =  s3X3(s)  [s2X(s)  -  sx(0)  -  x(0)]  (3-27) 

Comparison  of  equations  (3-19),   (3-25)  and  (3-27)  provides  the 
general  form  of  the  transform  for  these  functions. 


£[k(t)  xn(t)]  =   sn  Xn(s)  [s2X(s)  -  sx(0)  -  x(0)] 
=  snXn(s)   {X[x(t)]} 


(3-28) 
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C.  TRANSFORM  OF  THE  PRODUCT  OF  POWERS  OF  A  FUNCTION  AND  ITS 
DERIVATIVES 

Combining  the  results  of  the  two  foregoing  sections  the  general  form 
of  the  transform  of  the  product  of  a  function  raised  to  the  ntn   power  and 
its  various  derivatives  can  be  derived.    From  the  derived  transforms  as 
expressed  in  equations  (3-13)  and  (3-28): 

£[x(t)  Xn(t)]  =   snXn(s)   {£[x(t)]} 

X[x(t)xn(t)]  =   snXn(s)    U[x(t)]} 
The  Laplace  Transform  of  the  product  of  a  function  raised  to  the  ntn  power 
and  its  mtn  derivative  is  therefore: 

£[x(m)(t)xn(t)]  =  snXn(s)  (£[x(m)(t)]  }  (3-29) 

D,  THE  TRANSFORM  OF  THE  FIRST  DERIVATIVE  RAISED  TO  THE  nth  POWER 
From  Section  II-B,  equation  (2-5)  gives  the  transform  of  x(t). 
£[x(t)]  =  sX(s)  -  x(0) 

Proceeding  with  the  first  derivative  squared,   [x(t)]2,  and  expressing 
the  transform  using  integration  by  parts,  the  Laplace  Transform  becomes: 


03 


£  L(x)2]=  fx2  e~st=  x  fxe~stdt-x  []'xe"stdtdt  +  ..  .  (3-30) 

J  J  «J 

0  0  0 

Knowing  the  transform  of  x(t)  the  terms  of  the  series  may  be  expressed 
separately. 

00 

Ix   =  xjxe~stdt  =  x(0)[sX(s)  -  x(0)]  (3-31) 


Is  =   -x 


'Jxe-stdtdt  =^^  [sX(s)  -  x(0)]  (3-32) 


14 


j  xe_stdtdtdt  =^T  LsX(s)  -  x(0)] 


(3-33) 


£  {  [x(t)'J2  }  =  |=1  In  =  Ix   +  Ia  +  I3    +    .  . . 


=    [sX(s)-x(0)][x(0)  +^+^  +  ...] 

=    [sX(s)-x(0)]^-f  +f  +f  +f  + 


=  s2  [sX(s)  -  x(0)]    [x(s) 
=  s[sX(s)  -  x(0)]2 


X 


(0) 


(3-34) 


Using  similiar  arguments  the  transform  of  [x(t)]     may  be  derived. 


05 


£  (x(t)]3  }  =  Jx3e-Stdt  =  xjx3e'stdt  -k'Jj'x2e"st  dtdt  +  . . .  (3-35) 

0  0  e 

Knowing  the  transform  of  [x(t)]2  from  equation  (3-34)  let: 


Ix  =  x  jx2e~stdt  =  x(0)  s  [sX(s)  -  x(0)]: 


(3-36) 


U  =  -  x 


no 


rjx2e   stdtdt  =     ^  s  [sX(s)  -  x(0)]! 

*3 


(3-37) 


00 
I3=xfjTx2e-stdtdtdt  =^  s  [sX(s)  -x(0)]2 


(3-38) 


Summing  the  above  terms 


£{[x(t)]3}  =  £=1  Ix  +is  +  I3  +  ... 

=  s[sX(.s)  -x(0)]2[x(0)  +^  +^  +  ...  ] 
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=  s3[sX(s)  -  x(0)]~ 


xlO}+x(0)+xiO}  +  x(0]  + 


=  s*  [sX(s)  -  x(0)]3  (3-39) 

From  equations  (3-34)  and  (3-39)  the  transform  of  [x(t)]n  may  be  obtained 


£{[x(t)]n}=sn-1  [sX(s)  -x(0)]n 


=  s11"1 


Utx(t)]} 


n 


(3-40) 
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IV.    APPLICATION  OF  LAPLACE  TRANSFORMS  TO  NON-LINEAR  SYSTEMS 

A.      EXAMPLE  OF  A  FIRST-ORDER  NON-LINEAR  SYSTEM,  EQUATION  1 

The  first  non-linear  system  to  which  the  derived  Laplace  Transforms 
are  to  be  applied  is  one  which  has  previously  been  considered  by  Baycura 
[Ref.  3]  and  Brady  [Ref.  4]  in  the  initial  work  that  has  been  accomplished 
in  this  area.    The  differential  equation  to  be  considered  is  of  the  form: 

x  +  ax  +  bx2  =   0 

The  reasons  for  using  this  particular  equation,  aside  from  providing 
some  measure  of  continuity  in  previous  work,  are  that  it  represents  a 
stable  system  and  is  by  far  one  of  the  simplest  non-linear  equations 
with  which  to  work. 

In  this  initial  analysis  the  author  will  bring  to  light  several  concepts 
and  techniques  that  utilize  Laplace  Transforms  to  derive  solutions  to 
various  equations  and  then  compare  these  to  solutions  obtained  from  more 
"conventional"  techniques.    The  general  concept  will  be  to  express  the 
transform  of  the  unknown  variable  in  a  series  of  terms  in  s,  the  solution 
to  which  may  easily  be  obtained  on  the  digital  computer.    The  two  solu- 
tions used  for  comparison  will  have  been  obtained  from  an  analytical 
method  and  by  computer  using  Runge-Kutta  fourth-order  integration  methods 
using  Adams -Moulton  Predictor  Corrector  with  error  check. 

Utilizing,  therefore,  the  equation: 

x  +  ax  +  bx3  =  0  (4-1) 

and  applying  the  derived  Laplace  Transforms  we  obtain: 
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bsX2(s)  +  (s  +  a)X(s)  -x(0)  =  0  (4-2) 

Letting  x(0)  =  c  then: 

bsX2(s)  +  (s  +  a)X(s)  -  c  =  0  (4-3) 

Using  the  quadratic  formula,  solve  for  X(s). 


„ ,  .         s  +  a     f  r ,       4bc s      ~\  c  "1 

x(s)=  1ST  l-^E^l^)5]    I 

The  radical  may  be  expanded  in  a  binomial  series  of  the  form  of: 


(4-4) 


(l+^l+nM  +  ^^f^1     +    "(n-l)f-2),f  +  . , .  (4.5) 

Letting  a  =  2,  b  =  1,  and  x(0)  =  c  =  1  the  series  expression  for  X(s)  is: 

,  .    _     1  s  2ss  5s3  14s4    _     42s&  132s6  _ 

X(Sj  "  s+2        (s+2)3  +  (s+2)5  '  (s+2)7  +  (s+2)@        (s+2)u    +  (s+2)13  '  "  ' 

(4-6) 
To  obtain  x(t)  the  inverse  transform  of  the  individual  terms  may  be  obtain- 
ed and  summed.    An  alternate  and  by  far  a  more  convenient  and  simpler 
method  was,  however,  used  in  which  X(s)  was  expressed  as  a  ratio  of 
polynomials  in  s.    X(s)  then  took  the  form: 

.      n-1       .      n-2  , 

X(s)  =  bis  +  b2s  +  -  • ;  ■  +  bn 

sn  +  a^n-1  +  a2sn"z  +.  „  .  .  +  an 

A  devised  computer  technique  was  then  used  to  obtain  x(t).    This  method 

to  find  the  inverse  Laplace  Transform  by  digital  computer  is  thoroughly 

described  in  Appendix  A. 

Using  these  computer  methods  a  solution  to  the  equation  was  easily 

obtained  using  any  desired  number  of  terms  of  the  series.    As  a  result, 

solutions  to  the  equation  were  examined  terminating  with  the  second  to 

the  eighth  term  in  the  series. 
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Several  sets  of  solutions  have  therefore  been  derived  and  are  tabulat- 
ed in  Table  IV- 1.    These  are  an  exact  solution  which  is  known  and  obtain- 
ed through  analytical  methods  [Ref.  5],  a  computer  solution  using  Runga- 
Kutta  integration  methods,  and  several  computer  solutions  using  Laplace 
Transforms  but  using  varied  numbers  of  terms  in  the  truncated  series.     The 
Runga-Kutta  solution  is  presented  in  order  to  show  its  accuracy  since  this 
solution  will  be  the  basis  for  comparison  in  subsequent  examples. 


TABLE  IV- 

■1 

SOLUTION  TO  EQUATION  1 

EXACT 

RUNGE-KUTTA 

APPROX. 

APPROX. 

APPROX. 

TIME(t) 

SOLUTION 

SOLUTION 

4  TERMS 

6  TERMS 

8  TERMS 

0.0 

1.000 

1.000 

1.000 

1.000 

1.000 

0.2 

0.576 

0.575 

0.579 

0.579 

0.582 

0.4 

0.352 

0.352 

0.367 

0.369 

0.370 

0.6 

0.223 

0.223 

0.247 

0.249 

0.252 

1.0 

0.0950 

0.0944 

0.131 

0.101 

0.140 

1.4 

0.0410 

0.0414 

0.0841 

-0.0761 

0.0959 

1.8 

0.0184 

0.0183 

0.0574 

-0.464 

0.0791 

2.2 

0.00823 

0.00822 

0.0370 

-1.231 

0.0939 

2.6 

0.00368 

0.00368 

0.0216 

-2.434 

0.189 

3.0 

0.00165 

0.00165 

0.0116 

-3.969 

0.521 

3.4 

0.000734 

0.000735 

0.00636 

-5.594 

1.402 

3.8 

0*000333 

0.000327 

0.00416 

-7.024 

3.246 

4.2 

0.000153 

0.000145 

0.00350 

-8.028 

6.407 

4.6 

0.0000667 

0.0000643 

0.00334 

-8.482 

11.003 

5.0 

0.0000300 

0.0000283 

0.00318 

-8.387 

16.817 
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In  comparing  these  various  solutions,  it  is  noted  that  the  approximate 
transform  solution  conforms  quite  favorably  with  the  exact  solution  using 
up  through  five  terms  of  the  series  of  equation  (4-6) .    Truncating  beyond 
five  terms,  the  series  has  a  tendency  to  diverge,  the  direction  and  mag- 
nitude of  divergence  depending  on  the  number  of  terms  beyond  five  that 
are  utilized.    Since  this  series  should  converge  for  all  values  of  t,  it 
can  be  concluded  that  the  original  derivation  of  the  transform: 

£  [x(t)]  =  — ^    +       2~    +   ~3    +  ...  is  invalid  for  large  values  of  t  and 
is  the  direct  cause  of  this  divergence.    For  small  values  of  time,    t, 
however,  all  solutions  closely  approximate  the  true  solution. 

B.      EXAMPLE  OF  A  NON-LINEAR  EQUATION  WITH  SQUARE  OF  DERIVATIVE, 
EQUATION  2 

Considered  now  is  the  Froude  Equation  in  its  general  form: 

x  +  2£(l+ax)x  +  x=0  (4-7) 

This  equation  demonstrates  non-linear  damping,  but  otherwise  behaves 

in  a  rather  linear  fashion.    For  a  stable  system  select  £  =  0.1,  a  =  0  .5 . 

Then  equation  (4-7)  becomes: 

x   +    .2x  +    .  l(x)2  +  x  =   0  (4-8) 

Taking  transforms  we  have: 

s2X(s)  -  sx(0)  -  x(0)  +  .2sX(s)  -  ,2x(0)  +  .ls^C3(s)  -  .2x(0)s2X(s) 

+    .l[x(0)]2s  +X(s)  =   0.  (4-9) 

Collecting  terms  and  letting  x(0)  =  1  and  x(0)  =  0, 

.ls3Xs(s)  +  (.8s2  +  .2s  +  l]X(s)  -  (.9s  +  .2)  =  0  (4-10) 
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Now  solve  for  X(s)  using  the  quadratic  formula  and  arrange  in  a 
suitable  form: 

x(s).(.8S^.2s+1){,1±[l  +  ,SSs;.^]t}  (4.u, 

Expanding  the  term  within  the  radical  in  a  binomial  series  and  after 
collecting  terms  and  rearranging,  X(s)  can  be  expressed  in  a  series. 

-f  )  -         ('9s  +  *2)  .ls^(.9s  +  .2)s  ,02s6(.9s  +  .2)3 

MS)  "  (.8s2  +  .2s  +  1)        (.8sa+.2s  +  l)3        (.8s2+  .2s +  1)5 

(4-12) 
The  inverse  Laplace  Transform ,  x(t) ,  is  now  obtained  using  the 
method  depicted  in  Appendix  A.    The  solution,  compared  with  the  equa- 
tion solution  using  a  Runge-Kutta  integration  method,  is  displayed  in 
Table  !'.■'•  °  . 

The  Runge-Kutta  (true)  solution  and  the  approximate  solution  with 
values  corresponding  to  those  displayed  in  Table  IV-2  are  depicted  in 
graphical  form  in  Figure  1.    The  true  solution  is  a  damped  sinusoid.     The 
approximate  solution  in  this  case  is  fairly  accurate  for  small  values  of 
t  (up  to  3  seconds).    This  conforms  with  the  analysis  previously  given  in 
Section  IV-A. 

Introducing  appropriate  initial  conditions  at  periodic  intervals  of 
approximately  three  seconds,  several  piecewise  but  continuing  solutions 
were  obtained.    This  method  provided  a  continuous  solution  but  one  that 
could  be  obtained  over  short  periods  of  time.     Figure  1  shows  that  a  very 
accurate  solution  could  be  obtained  in  this  manner. 
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TABLE  IV-2 


SOLUTION  TO  EQUATION  2 


TIME 

0.0 

0.4 

0.8 

1.2 

1.6 

2.0 

2.4 

2.8 

3.2 

3.6 

4.0 

4.4 

4.8 

5.2 

5.6 

6.0 

6.4 

6.8 

7.2 

7.6 

8.0 


RUNGE-KUTTA 

TRANS      FORM 

SOLUTION 

SOLUTION 

1.000 

1.000 

0.922 

0.899 

0.709 

0.735 

0.397 

0.466 

0.038 

0.121 

-0.311 

-0.253 

-0.596 

-0.595 

-0.774 

-0.769 

-0.823 

-0.807 

-0.743 

-0.707 

-0.558 

-0.493 

-0.308 

-0.207 

-0.037 

0.054 

0.212 

0.205 

0.405 

0.400 

0.520 

0.520 

0.546 

0.557 

0.487 

0.515 

0.358 

0.407 

0.182 

0.249 

-0.012 

0.061 
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In  contrast  to  this  method  Figure  2  depicts  a  true  and  approximate 
solution,  the  latter  being  obtained  in  a  continuous  fashion  from  zero  to 
eight  seconds.     It  is  noted  that  prior  to  three  seconds  of  time  the  solution 
is  fairly  accurate  but  subsequent  to  this  time  begins  to  diverge  and 
becomes  slightly  out  of  phase.    The  divergence  is  noted  to  increase  with 
time . 

C.     EXAMPLE  OF  THE  PRODUCT  OF  THE  VARIABLE  AND  ITS  DERIVATIVE, 
EQUATION  3 

We  next  consider  the  van  der  Pol  Equation  which  describes  various 
physical  situations,  the  more  noted  representing  an  electronic  oscillator. 
The  general  expression  for  this  equation  is: 

x"  -  a  (1  -x2)x  +  kx  =  U 
Where  U  represents  the  input  and  a.  and  k  are  contants.    Consider  the 
system  with  zero  input  and  with  k=l,  then  the  equation  becomes: 

x"  -  a  x  +  a  x3x  +  x  =  0  (4-12) 

Apply  derived  transforms. 

s3  X(s)  -  sx(0)  -  x(0)  -as  X(s)  +  a  x(0)  +  a  s3X3(s)  -  a  s2x(0)Xa(s) 

+  X(s)   =   0  (4-13) 

After  collecting  terms  the  expression  becomes: 

as^sj-Qs^^s   +   [s3-as+l]X(s)-[sx(0)-ax(0)  +x(0)]=0 

(4-14) 

Solving  this  cubic  equation  for  X(s)  is  no  easy  task;  however,  if  we 

assume  that  the  initial  conditions  for  the  higher  derivative  terms  are 

negligible  and  equally  zero  then: 
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X(s)  =  2i°)    +    m  (4_15) 

s  s 

Substituting  for  X(s)  in  the  cubed  term  yields: 

asx(0)X2(s)  +  [s2  -  as  +  l]X(s)  -  [sx(0)  -  ax(0)  +  x(0)]  =  0         (4-16) 
Let:    x(0)  =  a 
x(0)  =  b 

x(0)  -  ax(0)  =  c 
Substituting  these  values  in  the  equation  for  X(s)  a  quadratic  equation 
is  obtained: 

absX2(s)  +  [s2  -  s  +  l]X(s)  -  [as  +  c]  =  0  .  (4-17) 

Using  the  quadratic  formula  now  solve  for  X(s). 

l 
X(s)   =  -(s2-QS    +   1)    +  [(s2  -  g  s    +   l)2  +   (4bg  s)  (as  +  c)  J 

2abs 


2 

2abs  l~x-L"L    '      (s2-as  +  l)3      JJ  (4-18) 


r  ffljp)    (as+c)    i| 

-   L  S2-ttS  +  l)2        JJ 


Expanding  the  term  under  the  radical  in  a  binomial  series  and  collecting 
terms  yields: 

yt  \        (as  +  c)  (bgs)(as  +  c)2        2(bgs)2(as  +  c)a 

K{S)  (S2-QS    +    1)  (S2_as+1)3+       {s3_as    +    1)5-  ... 

(4-19) 
Now  the  original  equation,   (4-11),  represents  a  system  having  variable 
damping  depending  on  the  displacement  x  and  hence  the  coefficient  of  the 
x  term.    Allowing: 

a  =  1 

x(0)  =  a  =  1 

x(0)  =  b  =  1 

x  -  ax(0)  =  c  =  0 
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equation  (4-19)  becomes: 

s  s3  2s5 

X(S)   =    (s2  -  s  +  1)       "    (s2  -  s  +  l)3     +    (s2  -  s  +  l)5     "  ••••     (4"20) 

Again  solve  for  x(t)  using  the  devised  computer  methods.    The  results  are 

compared  with  a  Runge-Kutta  solution  in  Table  IV-3.    For  a  small,  as  in 

our  example,  a  number  of  cycles  of  oscillation  are  required  before  the 

steady  state  is  achieved.    The  true  solution  is  shown  over  a  period  of  one 

cycle  of  the  sinusoidal  solution  but  during  this  period  the  initial  growth 

can  be  detected  commensurate  with  the  slight  increase  in  amplitude  with 

appropriate  increase  in  time. 

It  should  be  noted  at  this  time  that  the  derived  series  representation 

for  X(s)  as  depicted  by  equation  (4-20)  is  an  unstable  system.    This  fact 

is  borne  out  by  the  approximate  solution  as  shown  in  Table  IV-3.    It  can 

be  generally  stated  then  that  the  technique  adopted  here  for  finding 

Laplace  Transforms  is  not  applicable  to  systems,  such  as  represented  by 

the  van  der  Pol  Equation,  which  change  from  a  condition  of    instability  to 

one  of  stability.    The  method  being  investigated  will  not  show  the  true 

system  but  will  cause  the  solution  to  diverge  without  bound  with 

increasing  time . 
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TABLE  IV -3 
SOLUTION  TO  EQUATION  3 


TIME 
0.0 
0.2 
0.4 
0.6 
0.8 
1.2 
1.6 
2.0 
2.4 
2.8 
3.2 
3.6 
4.0 
4.4 
4.8 
5.2 
5.6 
6.0 
6.4 

6.8 
7.2 


RUNGE-KUTTA 

TRANS      FORM 

SOLUTION 

SOLUTION 

1.000 

1.000 

1.176 

1.174 

1.295 

1.273 

1.351 

1.256 

1.349 

1.079 

1.205 

0.071 

0.902 

-2.049 

0.421 

-5.374 

-0.319 

-9.583 

-1.261 

-13.723 

-1.850 

-16.041 

-1.913 

-13.938 

-1.744 

-4.184 

-1.473 

16.526 

-1.109 

50.561 

-0.601 

97.907 

0.160 

154.297 

1.187 

209.379 

1.893 

245.178 

1.994 

235.719 

1.837 

148.605 

TH 


V.     CONCLUSION 

In  the  foregoing  analysis  a  logical  method  was  presented  by  which 
the  Laplace  integral  was  expanded  into  an  infinite  series  and  after  inte- 
grating by  parts  the  infinite  series  developed  could  then  be  applied  to 
non-linear  systems.    It  was  shown  that  this  technique  is  applicable  to 
most  systems  or  differential  equations,  the  exceptions  being  systems  that 
are  unstable.    The  limitations  restricting  the  use  of  this  method  arise 
primarily  from  the  complexity  of  the  equation  and  the  inherent  difficulty 
in  solving  for  X(s),  the  transform  of  the  dependent  variable  x(t).    The 
procedure  used  in  this  paper  was  to  restrict  the  expression  for  X(s)  to 
a  quadratic  equation  and  to  further  express  the  transform  in  a  converging 
series  using  a  binomial  series  expansion. 

It  was  observed  that  in  obtaining  the  inverse  of  the  transform  a  close' 
approximation  to  the  true  equation/system  solution  could  be  obtained. 
This  approximation  was  found  to  be  valid  for  small  values  (several  sec- 
onds) of  solution  time  but  for  larger  values  of  time  a  slow  digression  from 
the  true  solution  was  obtained.    This  anomaly  was  attributed  to  the  errors 
introduced  through  the  use  of  the  transform  expressed  in  a  series  of  terms 
consisting  of  the  initial  conditions  of  the  unknown  variable  and  taking  the 
form  of: 


x(s)  =  *io)  +  *f  +  Mio)  +  _ 


In  the  attempt  to  find  a  reliable  yet  uncomplicated  method  to  produce 
the  inverse  of  the  Laplace  Transform  and  hence  the  solution  to  the  system, 
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a  computer  program  was  developed  whereby  these  painstaking  computa- 
tions could  then  be  readily  handled  on  the  digital  computer.    This  pro- 
cedure permitted  a  more  detailed  analysis  of  the  problem  otherwise 
unobtainable  through  ordinary  analytical  methods. 

Finally,  since  the  approximate  solution  was  found  to  be  valid  for 
only  a  small  interval  of  time  a  more  meaningful  comparison  with  the  true 
solution  was  necessary.    Through  the  introduction  of  appropriate  initial 
conditions  at  periodic  intervals,  several  piecewise  yet  continuing  solu- 
tions to  an  example  equation  were  obtained.    This  was  easily  executed 
with  the  use  of  the  associated  digital  program.    The  solution  obtained 
in  this  fashion  proved  very  accurate  and  usable  for  all  time. 
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APPENDIX  A 
COMPUTATION  OF  INVERSE  LAPLACE  TRANSFORM 

The  solution  of  example  non-linear  equations  through  application  of 
derived  Laplace  Transforms  is  based  primarily  on  computational  techniques 
proposed  by  Ward  and  Strum  [Ref.   6].    These  techniques  are  based  on 
expressing  the  system's  equations  in  simultaneous  first-order  state  equa- 
tions such  as: 

*(t)  =  Ax(t)  +  Bu(t)  (A-l) 

where:    x(t)  is  the  matrix  of  state  variables 

and  u(t)  is  the  matrix  of  system  inputs. 

Incorporating  a  time  domain  solution  of  these  state  equations,  x(t), 

computed  at  any  particular  time,  t,  may  be  then  expressed  as: 

t 

x(t)  =  c^1  x(0)  +  e^  [  €~At  By    (7)  dT  (A_2) 

■ 

0 

A  computer  program  is  presented  by  the  authors  to  solve  for  x. 
With  the  transform  of  x  given  as: 

,       n-1      ,       n-2 

X(s)  =  A! thl +—  -+V (A-3) 

sn  +  ax  sn"     +  a2sn       + +  an 

the  problem  then  is  to  transform  equation  (A-3)  into  the  form  of  equation 

(A-2). 

By  cross  multiplying 

,  n  n-1  n-2  .  „y  *      .       n-1      .       n-2  . 

(s     +  ax  s  +  a2s  +  ....  +  an)  X(s)  =  bx  s  +  b2s         +  .  .  .  .  +  bn 

(A-4) 
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The  right  side  is  then  written  in  matrix  form 


P  n-1        n-2        n-3  .  -, 

|_s  s  s  ....   1 J 


b2 


(A-5) 


If  X(s)  resulted  from  the  differential  equation: 


(n) 


x       +  a1  x 


(n-1)  (n-2)  .      ;■!.    .    ...    .  ._. 

+  a2s  + +anx=0  with  initial  conditions 


x(0)  ,x        (0) , x  (0)  then  the  transform  of  the  equation  is: 

snX(s)  -  sn_1x«»  -  sn"2x(l)  (0)  - -  s°x(n"1,(0) 


P  n-1  Tl  ,         n-2    ,_.         n-3    (i)   ,_. 
+  ax[s        X(s)  -  s        x(0)  -  s        xx      (0)  - 


-  s°x(n-2)  (0)] 


•    •••*• 


+  a      .  [sX(s)  -  s°x(0)] 
n-1 


+  a  X(s)  =  0 


(A- 6) 


Now  collecting  all  I.C.  terms  and  moving  them  to  the  right,  the 
right  side  of  the  equation  then  becomes: 

r  n-1  n-2  n-3  -,     ,nS 

=  |_s         +  &i  s         +  a2s         +...,+  a^^j  x(0) 

r  n-2              n-3  -,     (i)  ,„. 

+  [s  +a:s         + an_2  ]  xv    (0) 


+ 


+'G'iiei"1>(tt 


(A-7) 
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Rearranging  again  into  a  matrix  product: 


r  n-1      n-2 

=  [s  s 


s     1] 


1 


a      na      0.      .      .».      1 

n-1       n-2 


x(0) 
x(1)(0) 


x(n-1}(0) 


(A- 8) 


Now  equations   (A-5)  and  (A-8)  may  be  equated  so  that: 


ba 


n 


a  a 

n-1       n-2 


x(0) 

x(1)(0) 

x(2)(0) 


x'^V) 


(A-9) 


For  easier  notation  the  above  matrices  may  be  written  as: 

b  =  a  x(0)  (A- 10) 

or:    x(0)  =a_1b  (A- 11) 

Applying  these  results  to  equation  (A-2)  note  that  in  this  case  the 

input,  the  u  matrix,  will  be  zero,  leaving  but  one  term  for  x(t). 

x(t)  =  e~l  x(0)  (A-12) 

where  A  is  a  matrix  of  constants  as  given  in  equation  (A-l)  and  is  of  the 
form: 
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-a 


-a 


-a 


n  n-1 1 

Knowing  a  in  equation  (A-ll)   x(0)  may  be  easily  obtained  with  a  compu- 
ter  utilizing  a  library  program  to  obtain  the  inverse  of  a  matrix  and  then 
multiplying  the  two  resulting  matrices. 
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APPENDIX  B 
LIST  OF  LAPLACE  TRANSFORMS 

1.  £[x(t)]  =  X(s) 

2.  X[x2(t)]  =   sX2(s) 

3.  X[x3(t)]  =  s*X3(s) 

4.  I[xn(t)]=  sn"1Xn(s) 

5.  JE[x(t)]  =   sX(s)  -  x(0) 

6.  £[x"(t)]  =   s2X(s)  -  s  x(0)  -  x(0) 

7.  X[xv    y(t)]  =   s    X(s)    -   s  x(0)  -  s  x(0)  -  ....  -  sxv         '(0) 

8.  X[x(t)]  =  sX(s)    -  x(0) 

9.  £[{*(t)}2]  =  s[sX(s)  -x(0)]2 

10.  £[{x(t)}n]=  sn"1[sX(s)  -x(0)]n  =   sn'1    [X[x(t)}n 

11.  Zi^h)}*^**-1    [X[x(m)(t)]}n 

12.  x[x(t)  x(t)]  =   sX(s)   [sX(s)  -  x(0)] 

13.  £[x(t)  x2(t)]  =   s2X2(s)  [sX(s)  -  x(0)] 

14.  £[x(t)  xn(t)]  =  snXn(s)  U[x(t)]} 

15.  £[x*(t)  x(t)]  =  sX(s)  [saX(s)  -  sx(0)  -  x(0)] 

16.  JC[x(t)  x2(t)]  =   s2X2  s[s2X(s)  -  s  x(0)  -  x(0)] 

17.  £[x(t)xn(t)]  =   snxn(s)   {X[x(t)]} 

18.  X[x(m)(t)xn(t)]=   snXn(s)    {X[x(m)(t)]} 
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